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embryos 
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Abstract 

Background: Previous studies have suggested that the axon guidance proteins Slitl and Slit2 co-operate to establish 
the optic chiasm in its correct position at the ventral diencephalic midline. This is based on the observation that, 
although both Slitl and Slit2 are expressed around the ventral midline, mice defective in either gene alone exhibit few 
or no axon guidance defects at the optic chiasm whereas embryos lacking both Slitl and Slit2 develop a large 
additional chiasm anterior to the chiasm's normal position. Here we used steerable-filters to quantify key properties of 
the population of axons at the chiasm in wild-type, S//f7 -/ ~ S//f2~ /_ and S//f7 _/_ 5//f2 _/_ embryos. 

Results: We applied the steerable-filter algorithm successfully to images of embryonic retinal axons labelled from a 
single eye shortly after they have crossed the midline. We combined data from multiple embryos of the same 
genotype and made statistical comparisons of axonal distributions, orientations and curvatures between genotype 
groups. We compared data from the analysis of axons with data on the expression of Slitl and Slit2. The results showed 
a misorientation and a corresponding anterior shift in the position of many axons at the chiasm of both S//'f2 _/_ and 
S//'f 7 _/_ S//'f2 _/_ mutants. There were very few axon defects at the chiasm of S//t7 _/ " mutants. 

Conclusions: We found defects of the chiasms of S//'f7" / "S//'f2" / " and Slit 7 _/ " mutants similar to those reported 
previously. In addition, we discovered previously unreported defects resulting from loss of Slit2 alone. This indicates the 
value of a quantitative approach to complex pathway analysis and shows that Slit2 can act alone to control aspects of 
retinal axon routing across the ventral diencephalic midline. 
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Background 

The development of the complex connectivity of the ner- 
vous system involves the concomitant elongation and 
guidance of axons to specific targets. Growing axons are 
guided at their distal tips by growth cones, which move 
using relatively autonomous mechanisms to sense and re- 
spond to environmental cues. Some of these cues are 
molecules distributed in gradients that either attract spe- 
cific growth cones up their concentration gradients (che- 
moattractants) or repel growth cones away from their 
sites of increasing concentrations (chemorepellents). One 



* Correspondence: David.Price@ed.ac.uk 

2 Centre for Integrative Physiology, University of Edinburgh, Edinburgh, UK 
Full list of author information is available at the end of the article 



structure whose formation has provided an excellent 
model in which to study axon guidance is the optic 
chiasm of the mammalian visual system. 

Retinal ganglion cell (RGC) axons from the two eyes 
converge on each other at the midline of the ventral part 
of the brain (ventral to the hypothalamus) where they ei- 
ther cross the midline (the contralateral tract) or turn 
away from it (the ipsilateral tract), forming the x -shaped 
optic chiasm. The formation of this pathway is achieved 
through the growth of RGC axons that exit the retina in 
tight bundles to reach the midline at a position ventral 
to the hypothalamus, where they either cross or do not 
cross depending on the retinal locations of their cell 
bodies. Guidance molecules that have been implicated in 
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these regulatory tasks include Shh, which acts as a mid- 
line repellent to RGC axons except at the point where 
the chiasm forms [1], EphrinB2 and its receptor EphBl, 
which are critical for ipsilateral/contralateral sorting [2], 
and the extracellular secreted molecules Slitl and Slit2, 
which signal via Robo receptors to constrain axons to 
the region of the chiasm [3-6] . 

An interesting aspect of Slitl and Slit2 function at the 
chiasm is that, although both are expressed at the hypo- 
thalamic ventral midline and RGC growth cones respond 
to individual Slits [7-9], mice defective in either gene alone 
have been described as exhibiting either few or no RGC 
axon guidance defects at the optic chiasm [10]. By con- 
trast, in double-mutant mice lacking both Slitl and Slit2, 
a large additional chiasm develops anterior to the chiasms 
normal position [10]. These results indicate that Slit pro- 
teins co-operate to establish a corridor through which the 
axons are channelled, thereby helping define the site in 
the ventral diencephalon where the optic chiasm forms. 

The nature of the co-operation between Slitl and Slit2 
at the optic chiasm, for example whether one plays a 
stronger a role than the other, is not clear. One of the 
main reasons is that previous analyses of the structure of 
the optic chiasm and the routes taken by the axons it 
contains rely on non-quantitative anatomical descrip- 
tions of its appearances after most or all of its axons 
have been labelled. The complexity of the chiasm and 



the large number of axons that it contains make it diffi- 
cult to carry out quantifications that might reveal effects 
of the Slits individually. In this work we adapted a 
method based on second-derivative Gaussian steerable 
filters [11] to allow us to obtain quantitative information 
on the positions, curvatures and orientations of develop- 
ing axons in large tracts such as the optic chiasm. Steer- 
able filters have been used successfully to quantify key 
properties of neuronal processes in vitro [12], where the 
separation of the processes makes them clearer and eas- 
ier to analyze than in vivo. Using a steerable-filters- 
based technique tailored specifically for our purposes, 
we found previously unreported defects at the optic 
chiasm of Slit2~ x ~ single-mutant mice that correlated 
with the normal expression pattern of Slit2. Our results 
suggest that Slit2 makes a greater contribution than Slitl 
to the guidance of RGC axons at the optic chiasm. The 
application of a similar approach to the analysis of other 
mutants and similarly complex pathways might be a 
profitable way of finding hitherto undetected defects in 
some strains. 

Results 

Quantitative analysis of labelled axons at the wild-type 
chiasm 

We first examined Dil-labelled axons at the optic chiasm 
of wild-type mice (Figure 1). Figure 1A shows an example 





Figure 1 Applying steerable filters to a Dil image of the developing mouse optic chiasm. (A) A 150 um-thick horizontal section of an El 3.5 
mouse brain into whose left eye Dil had been injected. An image of the chiasm taken with fluorescence was superimposed, revealing the Dil 
labelled axons. A grid whose baseline ran between the centers of each eye was used to align images from multiple embryos. The red shaded 
rectangle highlights the region that was selected for statistical comparison in sets of mice of different genotypes. Scale bar = 200 urn. 
(B) Magnification of the Dil labelled chiasm from A. Scale bar = 200 urn. (C) 2nd-derivative Gaussian filters of different orientations used to 
analyze images: among these examples, the top filter would respond best along the axon in D. Scale bar = 1 urn. (D) High magnification of a 
segment of axon from B. Arrows, which are associated with each pixel along the centre of the axon and point away from the injected eye, are 
aligned with the filter orientations that gave the best responses. Scale bar = 2 urn. (E) The result of applying the full algorithm to the image in B: 
red arrows represent direction, green arrows represent curvature. Only a small subset of arrows has been plotted, for clarity. Scale bar = 200 urn. 
(F) Examples of two arrows at different positions on a single axon (surrounding axons removed for clarity) showing how vector length is 
proportional to curvature, with the direction of the arrow indicating the direction of curvature. Scale bar = 10 urn. 
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of the optic chiasm in a horizontal section at embryonic 
day 13.5 (E13.5). E13.5 was selected because retinal axons 
have recently crossed the midline at this age, having first 
penetrated the diencephalon at E12.5 Erskine et aL, [8]. It 
is, therefore, the earliest age offering the opportunity to 
observe potential defects of guidance at the chiasm in 
mutants. 

Figure IB is a high magnification view of the same 
chiasm as in Figure 1A, showing detail of the Dil la- 
belled axons. Axons such as these were analyzed using 
Gaussian steerable filters, examples of which are illu- 
strated as oriented spatial functions in Figure 1C. The 
filters were convolved with the images of Dil-labelled 
chiasms. Only data derived from points along axons 
were retained; points that were not on axons were 
excluded using non- maximum suppression, which 
removed data from areas that did not lie on ridges of 
high intensity with respect to the surrounding landscape 
of the image. Figure ID illustrates the outcome: the 
orientations of filter that gave the best response at each 
position along the axon are shown as arrows pointing 
away from the Dil-injected eye (referred to as direc- 
tions). An example of applying the algorithm to all the 
axons at the chiasm is illustrated in Figure IE, which 
shows only a small subset of directions (or oriented 
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vectors) for purposes of clarity. The subset is <0.01 of 
the size of the full set, allowing the arrows to be shown 
enlarged compared to those in Figure ID. The algorithm 
provides a method to quantify automatically a vector 
field representing the orientations of axons in a given 
image and also the curls of the vector fields, which 
measure the degree of curvature of axons at each point. 
Notice that the arrows representing direction (red 
arrows in Figure 1D,E) have a constant magnitude 
whereas those representing curvature (green arrows in 
Figure 1E,F) have both direction and length, with the 
length proportional to the amount of turning. 

Axon locations in Slit1 _/ ~ Slit2 _/ " and Slit1 _/ "; Slit2 _/ " embryos 

Figure 2 illustrates the results of applying the steerable 
filter algorithm to sets of E13.5 embryos that were wild- 
type (Figure 2A,E,I; n = 9 embryos), Slitl~'~ (Figure 2B, 
FJ; n = 11 embryos), Slit2~ f ~ (Figure 2C,G,K; n = 5 
embryos) or Slitl~'~Slit2~'~ (Figure 2D,H,L; n = 4 
embryos). Figure 2A-D shows individual examples of 
chiasms analyzed as outlined in Figure 1. As can be seen, 
there is considerable variation in the appearance of these 
examples, particularly between those that are Slit2~'~ or 
Slitr'~Slit2~'~ and those that are wild- type or Slitl'^. 
To assess whether this variation can be accounted for by 
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Figure 2 Combining data from multiple embryos with different Slit! and Slit2 genotypes. (A-D) Images depict examples of chiasms 
analyzed using the steerable filter algorithm, with conventions as in Figure 1E. Scale bars = 200 urn. (E-H) Plots show the direction vector fields 
overlaid (using the grid system in Figure 1A) from all embryos of each genotype, with each colour representing data from a different embryo. 
Only small subsets of arrows from each embryo have been plotted, for clarity. (I-L) Heat-maps: In each square we calculated the mean number of 
arrows, including all arrows from all embryos of each genotype (these values for mean densities of arrows are proportional to the mean densities 
of axon). A colour scale is used to illustrate the spatial distribution of mean densities on a relative scale for each genotype. A, E, I, wild-type, n = 9; 
B, F, J, S/ztr 7 -, n = 1 1; C,G,K, S//t2" / ", n = 5; D,H,L, Slitr'Slitr'', n = 4. 
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differences in genotype, we combined results from ap- 
plying the algorithm to all embryos of each genotype. 
Figure 2E-H shows examples of vectors (i.e. subsets, as 
described above) representing axonal directions from all 
embryos of each genotype. Data from different embryos 
are shown in different colours, providing a sample of the 
routes taken by axons in all embryos of each genotype. 
Data from embryos within each genotype group were 
aligned using the grid system shown in Figure 1A. The 
mean distance between the centre of the two eyes did 
not vary significantly with genotype, varying by < 8% be- 
tween groups. The plots in Figure 2E-H suggest an an- 
terior shift in the position of the optic chiasm in Slit2~'~ 
and Slitr'~Slit2~'~ embryos compared to wild- type and 
Slitr^ embryos. 

To test this statistically, the area of each chiasm was 
split into a 32 x 32 grid. Figure 2I-L shows the spatial 
distributions of the mean numbers of oriented vector 
fields (i.e. the total numbers of arrows, as exemplified in 
red in Figure ID) across the chiasms of all embryos of 
each genotypes (higher numbers are towards the red end 
of the spectrum). All vectors were included (i.e. not the 
samples explained above and used for illustrative pur- 
poses in Figures 1 and 2). These values were propor- 
tional to the mean densities of axon within each square 
for each genotype, since in all cases analysis with steer- 
able filters was done at constant intervals and filtering 
using non-maximum suppression prevented the inclu- 
sion of data from areas that contained no axons. Note 



that while these values are proportional to the densities 
of axon within each square, and can therefore be used to 
examine distributions of axon, they can not be used to 
derive values for the absolute numbers of individual 
axons across the chiasm. Our approach does not attempt 
to trace individual axons and can not, therefore, give 
their absolute numbers. 

The graphs in Figure 2I-L suggest an anteriorization 
of the population of chiasmatic axons in Slit2~'~ and 
Slitr l ~Slit2~'~ embryos. For each area of the chiasm we 
tested for significant differences (taking account of mul- 
tiple testing, see Methods) between the values from 
wild-type and Slitr 1 ', Slit2~'~ or Slitr'~Slit2~'~ 
embryos: results are plotted in Figure 3A,E,I. There were 
no significant differences between wild- type and Slitl~'~ 
embryos (Figure 3A), but significantly larger densities of 
axon were located in abnormally anterior positions in 
both Slit2~'~ and Slitr'~Slit2~'~ embryos (Figure 3E,I). 
There were no significant differences between Slit2~'~ 
and Slitr'~Slit2~'~ embryos (not shown). 

Axon orientations and curvatures in Slit1 _/ ~, Slit2 _/ ~ and 
Slitl" 7 "; Slit2 _/ " embryos 

The same procedure was then applied to values of orien- 
tation, curvature angle and curvature magnitude from 
each area across the chiasm. These values were the 
averages of the vectors within each square of the 32 x 
32 grid and, therefore, took no account of the trajector- 
ies of individual axons, e.g. whether they crossed each 
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Figure 3 Areas of statistically significant differences between wild-type and mutant chiasms. In each panel the white squares indicate 
positions that differ significantly from wild-type; Student's t-test was applied for axon number and Watson-Williams test for orientation and 
curvature. All tests were controlled for false positives [14]; a = 0.05; see Materials and Methods). Columns show differences in mean axon number, 
axon orientation, angle of axon curvature and magnitude of axon curvature per grid square. (A-D) S//'t7 _/ " embryos; (E-H) S//f2 _/ ~ embryos; 
(l-L) S//f7" / "S//f2" / " embryos. Wild-type, n = 9; 5//t7 _/ ~ n = 1 1; S\\X2~'~, n = 5; S\\xr'~S\\X2~'~ ", n = 4. 
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other or not. Comparison of wild- type and Slitl~ 
embryos showed only a few areas returning statistically 
significant differences in either axon orientation or angle 
of curvature (Figure 3B,C). Comparison of wild- type and 
either Slit2~ f ~ or Slitl' 1 ~Slit2~'~ embryos showed many 
more areas returning statistically significant differences 
in axon orientation, many of which were located lateral 
to the midline (Figure 3FJ; the midline runs vertically 
through the centre of each panel). In Slitr'~Slit2~ l ~ 
double-mutant embryos there were also significant ab- 
normalities of axonal orientations in a posterior area 
contralateral to the injected eye (areas in the bottom 
right of Figure 3J) that were not present in Slit2~'~ 
mutants. Regarding the angle of axon curvature, most 
differences between genotypes were found around the 
midline (Figure 3G,K). No comparison returned any sig- 
nificant differences in magnitudes of curvature 
(Figure 3D,H,L). Overall, these data indicate that, in 
Slit2~'~ and Slitr'~Slit2~ l ~ embryos, many axons are 
oriented abnormally in their route across the chiasm, an 
observation that agrees with our finding described above 
that many axons are mislocated. In Slitl~ f ~ embryos 
there were few axon orientation defects, in line with 
there being no detectable axonal mislocation. 

Comparison of axonal defects with patterns of Slit 
expression 

In agreement with previous studies, in situ hybridiza- 
tions at E13.5 revealed strong Slitl mRNA expression 
both anterior and posterior to the junction of the optic 
nerve and the brain (Figure 4A) [8] whereas expression 




Slit2 



Figure 4 Expression of Slitl and Slit2 around the El 3.5 optic 
chiasm. In situ hybridizations on horizontal sections through the 
ventral diencephalon: areas of expression appear pale yellow. 
Expression of Slitl is strong around the point of entry of the optic 
nerve (o.n.); expression of Slit2 is strong anterior to this point. 
Olfactory epithelium = o.e. Scale bar = 200 urn. 

V / 



of Slit2 was strongest anterior to the point of entry of ret- 
inal axons (Figure 4B) [7,8]. The analysis above revealed 
that many retinal axons of Slit2-'- and Slitr'SlitT 1 - 
mutants crossed at abnormally anterior locations and here 
we examined the spatial relationship between the orienta- 
tions of these axons and the normal Slit2 expression pattern. 

Using the same strategy described above for analysis of 
axons, we combined data from three separate comparably- 
developed in situ hybridizations to give a map of the aver- 
age staining intensity for Slit2 across the ventral midline in 
150 (im horizontal sections at the level of retinal axonal 
entry. The system illustrated in Figure 1A was used to align 
data on gene expression in wild-types with data on retinal 
axons from wild-type or mutant embryos. Results are 
shown in Figure 5. The following vector-fields were 
obtained. For axons, mean orientations (red/magenta in 
Figure 5B-D) and mean directions and magnitudes of 
curvature (green in Figure 5B-D) were obtained within 
each square in the 32 x 32 grid (described above), using 
data from all embryos of each genotype. Also within each 
square, the mean vector representing the gradient of Slit2 
expression was calculated (yellow in Figure 5B-D; the 
lengths of the lines represent the magnitudes of the gradi- 
ent with arrowheads pointing from high to low intensity of 
label). The location of these vector fields in the brain is 
shown in Figure 5A. 

Axons were labelled from the left eye. In wild-types 
(Figure 5B) labelled axons about 200um to the left of the 
midline were oriented roughly 45-75° relative to the mid- 
line. They turned to run roughly orthogonal to the midline 
as they as they approached it (note the green vectors 
representing curvature concentrated near to the midline 
in Figure 5B). In taking this course, the axons were 
oriented roughly (± about 30°) orthogonal to the vectors 
representing the Slit2 gradient and, therefore, avoided the 
anterior region of high Slit2 expression (Figure 5B). In 
contrast, axons on the left of the midline in S7/£2 _/ ~ 
mutants were oriented roughly orthogonal to the midline 
throughout their approach, thereby entering the anterior 
areas where Slit2 would normally be expressed. Many 
turned at the midline to exit through the contralateral 
area of high Slit2 expression (Figure 5C). A similar pattern 
was observed in Slitl~ l ~Slit2~ l ~ mutants (Figure 5D). 

All these analyses together provide a consistent picture 
in which loss of either Slit2 alone or Slitl and Slit2 to- 
gether result in many retinal axons being misoriented on 
approach and exit from the ventral midline and many 
being located abnormally anteriorly, where Slit2 would 
normally be expressed. It appears that Slit2 is required 
to prevent retinal axons from taking this anterior route. 

Discussion 

In this study we have shown it is possible to apply rela- 
tively simple image analysis methods to static images of 
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Figure 5 Axonal defects in S//f2 _/ " and Slitl^'Slte''' embryos 
occur where Slit2 is normally expressed. (A) The area of axonal 
analysis is superimposed on an in situ hybridization for Slit2. Scale 
bar = 200 urn. (B-D) Mean vector orientation (red/magenta), 
direction and magnitude of curvature (green) and magnitude and 
direction of Slit2 gradient (yellow) are shown for each grid square for 
each genotype shown. Wild-type, n = 9; S//'t2" / ", n = 5; S//'t7" / "S//'t2" / " 
n = 4. Data are superimposed on an image showing the average 
intensities within each grid square of Slit2 expression in wild-types 
(n = 3 in situ hybridizations). 



developing axons, quantify key properties relevant to the 
study of axon guidance and combine data from multiple 
embryos to make statistical comparisons between groups 
sharing a particular property (here their genotype). Since 
it is very difficult to trace individual axons in dense 
axonal tracts in vivo, our approach was to determine the 
average positions, directions and curvatures of popula- 
tions of axons. Data were combined from multiple 
embryos so as to allow statistical comparisons that took 
into account variations between individual embryos or 
introduced by technical factors such as the degree of fix- 
ation, the amount of tracer injected or the precise plane 



of section in each individual. Using this approach, we 
detected the previously-described defect of chiasm de- 
velopment in Slitr'~Slit2~ l ~ mutants [10]. The real 
value of additional objectivity was demonstrated by our 
detection of hitherto unreported defects in Slit2~ f ~ 
mutants. 

Based on their qualitative examination of mutant 
chiasms, Plump et al. [10] stated that they were unable 
to detect defects at the chiasm of Slit2~'~ mutants, with 
the caveat that they could not exclude defects that were 
beyond the sensitivity of their experiments. In fact, in 
the example that they show of labelled axons at the 
chiasm of Slit2~'~ mutants (Figure 3 of their paper) the 
anteroposterior width of the tract at the midline is about 
double the width of the wild-type. Our quantitative results 
indicate that there is indeed a defect of the Slit2~'~ chiasm 
with many axons mislocated abnormally anteriorly. Our 
results agree with the conclusion of Plump et al. [10] that 
there are no major defects at the chiasm of Slitl' 1 ~ 
mutants. Regarding the comparison of Slit2~'~ and 
Slitr'~Slit2~'~ embryos, we found significant misor- 
ientation of some of the posterior axons contralateral 
to the injected eye in double mutants, but otherwise 
they were similar to single mutants. Such misorienta- 
tion of posterior axons was not found in Slitr'~ 
mutants, indicating that there are abnormalities of 
double mutants that are not present in either mutant 
alone. Overall, however, our results indicate that, early 
in the formation of the chiasm, Slit2 plays a more 
powerful role than Slitl in constraining the growth of 
axons to their correct location across the ventral 
midline. 

While our results suggest a greater similarity between 
the effects of loss of Slit2 alone and loss of both Slitl 
and Slit2 than was suggested by Plump et al. [10], the 
results of the two studies are in fact not strictly compar- 
able. Plump et al. [10] reported a major difference be- 
tween the Slit2~'~ and Slitr'~Slit2~'~ genotypes at 
E15.5, which is two days later than our findings. Al- 
though they reported defects of the double mutants at 
earlier ages, including E13.5, that look very similar to 
those found here, they did not provide data on Sy/£2 _/ ~ 
embryos at earlier ages. This raises the interesting possi- 
bility that the effects of losing both Slitl and Slit2 be- 
come progressively more severe than those of losing 
Slit2 alone as the chiasm develops from E13.5 to E15.5. 

The Slit'^ phenotypes reported here have a striking 
correlation to the ectopic projection of Slit~ x ~ retinal 
axons to the contralateral eye that we reported previ- 
ously. In wild- type and Slitl~ x ~ embryos the inter- 
retinal projection is relatively small indicating that 
Slitl is dispensable for keeping retinal axons out of the 
opposite eye. In contrast, loss of Slit2 function results 
in a dramatic increase in the size of the inter-retinal 
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misprojection and a further increase occurs in Slitl' ~; 
Slit2~ / ~ embryos [6]. It is easy to imagine that these phe- 
notypes are causally linked: those retinal axons that cross 
the midline in aberrant positions in Slit mutants are liable 
to end up on a track which predisposes them to mispro- 
ject to the opposite eye. This provides an example of how 
the novel analysis techniques described here can shed 
extra light on axon guidance phenotypes previously 
observed using more traditional techniques. 

The effect of Slit2 is likely explained by its spatial pat- 
tern of expression, which was studied here by in situ 
hybridization due to the lack of suitable antibodies. 
Clearly, it would be preferable to examine the distribu- 
tion of Slit proteins, and our analysis of this issue is 
based on the assumption that the protein distribution 
would approximate to the mRNA distribution at a tissue 
level. Slit2 mRNA is expressed in the ventral midline in 
a position anterior to the normal chiasm. It is straight- 
forward to understand how loss of this expression might 
be a critical factor allowing axons to cross in an abnor- 
mally anterior position. The effects of Slitl loss are intri- 
guing and less easily explained. Slitl is normally 
expressed both anterior and posterior to the point of 
entry of retinal axons. Its loss in combination with that 
of Slit2 causes some misorientation of posterior axons 
after they have crossed the midline, but why this defect 
is not detected in single Slitr x ~ embryos is not clear. It 
appears that the presence of other factors provides suffi- 
cient guidance even in the absence of Slitl and that Slit2 
might be one of these other factors. How Slitl and Slit2 
cooperate to prevent contralateral posterior misguidance 
is not clear. It is possible that Slit2 prevents axons 
approaching the midline from acquiring abnormalities 
that predispose some of them to require Slitl repulsion 
from posterior contralateral territory. In this scenario 
Slitl repulsion would only be required if Slit2 is lost. 

Conclusion 

Our results identify a previously undetected but import- 
ant role of Slit2 alone at the chiasm, while also showing 
some degree of cooperation between Slitl and Slit2 and 
a lack of an obvious role for Slitl alone, as suggested by 
previous work (Plump et al, 2002). We suggest that the 
approach adopted here can increase the sensitivity for 
detecting axonal defects in mutant strains and might be 
adopted more widely in the future. 

Methods 

Mice and mutant alleles 

The licence authorising this work was approved by the 
University of Edinburgh Ethical Review Committee of 
22 nd September 2008 (application number PL35-08) and 
by the Home Office on 6 th November 2008 (licence num- 
ber 60/3913). Animal husbandry was in accordance with 



the UK Animals (Scientific Procedures) Act 1986 regula- 
tions. Slitl and Slit2 mutant alleles were described by 
Plump et al. [10]. The Slitl' allele was identified by multi- 
plex PCR using primers 5-ACCCTTAGCTTCTACCAA 
CC-3, 5- TCTCCTTTGATCTGAGACCG-3 and 5-AGG 
TTTCTCGAGCGTCATAG-3: the wild-type allele gives a 
544bp product, the mutant allele gives a 393bp product. 
The Slit2' allele was identified by multiplex PCR using pri- 
mers 5-AAGACCTGTCGCTTCTGTCAG- 3, 5-AAACA 
GGTTTCTACCGCACG-3, and 5-AAGTCTAGTAGAGT 
CGAGCG-3: the wild-type allele gives a 600bp product, 
the mutant allele gives a 350bp product. All mice were 
C57BL6. 

Dil tract tracing 

E13.5 embryonic heads were fixed at 4°C in 4% parafor- 
maldehyde in phosphate buffered saline (PBS) overnight, 
and 1 , 1 -dioctadecyl-3,3,3,3-tetramethyl-indocarbocyanine 
perchlorate (Dil) crystals (Invitrogen, San Diego, CA) were 
placed into the optic cup of one eye, after removal of the 
lens, to label axons leaving the retina. Heads were 
returned to 4%paraformaldehyde in the dark at room 
temperature for 4 weeks to allow tracers to diffuse along 
axons. Heads were then sectioned (150 um) with a vibra- 
tome, cleared in 9:1 glycerol: PBS and mounted in Vecta- 
shield (Vector Laboratories, Burlingame, CA). Images 
were acquired at x40 using a Zeiss Axiovert confocal LSM 
510 microscope in tile-scan mode to collect serial optical 
sections, which were then combined to create a projec- 
tion. The section thickness was sufficient to include the 
entire chiasm but where, due to the position of the cuts, 
part of it was in a second section, both sections were 
imaged and a composite was generated so that the entire 
depth of the chiasm was analyzed in all cases. 

Image analysis 

Images were imported into MATLAB (Math Works, 
Natick, MA) and convolved with a bank of three 2D 
2nd-derivative Gaussian filters of width 1.8 um. The 
maximum filter response was calculated. The method is 
described in full in Freeman and Adelson [11]. Imaging 
was done with consistent confocal settings chosen to 
avoid saturation at high magnification. Rather than ad- 
just settings, where small areas of saturation did occur 
and the filters did not activate, vector orientations in the 
saturated area were interpolated based on the orienta- 
tions the vectors were taking on the edges of that area. 
In practice, this approach only applied in very small 
areas in some images. Non-maximum suppression was 
then performed by considering the two orthogonal pixels 
at any given point in the direction vector field and keep- 
ing it if both these flanking pixels gave a lower filter re- 
sponse. A global threshold was then applied that 
removed all pixels for which the intensity was less than 
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twice the median intensity across all pixels or whose fil- 
ter response was less than twice the median. The orien- 
tation vector field was then smoothed with a Gaussian 
of 30 \im to remove noise (which is of a higher fre- 
quency) and the curl of this vector field was calculated 
to give the curvature of the axon trajectories at any 
given point. Non-maximum suppression was applied a 
second time. 

Statistical comparison of genotypes 

To compare different genotypes with respect to axon 
number, axon direction and axon curvature across the 
image, the vector fields were rotated and scaled to a uni- 
versal grid with the left eye of the embryo defining the 
origin and the right eye defining the position (0, 1), as 
shown in Figure 1A. The area of the chiasm was split 
into a 32 x 32 grid covering the region (0.4 : 0.6, 0.05 : 
0.25), highlighted in Figure 1A. The axon number, mean 
axon orientation, mean curvature direction and mean 
curvature magnitude were put into each of these bins 
for every sample. A Student s £-test was applied for axon 
number and the Watson- Williams test was used to com- 
pare axon orientations and angle of curvature between 
samples, as provided by the CircStat MATLAB toolbox 
(circular statistics is appropriate for analysis of curva- 
tures and orientations) [13]. To control for false posi- 
tives, FDR (False Discovery Rate) [14] was used, given 
the large number of hypotheses present. For all tests we 
took a = 0.05. These methods are commonly used for 
dealing with multiple comparisons in other areas of 
spatial statistics. 

In situ hybridization 

Embryonic heads were fixed overnight in 4% paraformal- 
dehyde in phosphate buffered saline (pH 9.5 at 4°C). In 
situ hybridizations for Slitl or Slit2 used 100 \im vibra- 
tome sections and digoxigenin-labeled antisense ribop- 
robes synthesed from rat cDNAs encoding Slitl and Slit2 
as templates, as previously described Erskine et al, [8]. 
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